The Brown Dwarf Kinematics Project (BDKP) I. Proper Motions 
and Tangential Velocities for a Large Sample of Late-type M, L 

and T Dwarfs 

Jacqueline K. Faherty^ Adam J. Burgasscr^, Kclle L. Cniz'^'^, Michael M. Shara^, 
Frederick M. Walter^, and Christopher R. Gelino^ 

ABSTRACT 

We report proper motion measurements for 427 late-type M, L and T dwarfs, 
332 of which have been measured for the first time. Combining these new proper 
motions with previously published measurements yields a sample of 841 M7-T8 
dwarfs. We combined parallax measurements or calculated spectrophotometric 
distances and computed tangential velocities for the entire sample. We find that 
kinematics for the full and volume-limited 20 pc samples are consistent with 
those expected for the Galactic thin disk, with no significant differences between 
late-type M, L, and T dwarfs. Applying an age-velocity relation we conclude 
that the average kinematic age of the 20 pc sample of ultracool dwarfs is older 
than recent kinematic estimates and more consistent with age results calculated 
with population synthesis models. There is a statistically distinct population of 
high tangential velocity sources {Vtan > 100 km s^-*^) whose kinematics suggest 
an even older population of ultracool dwarfs belonging to cither the Galactic 
thick disk or halo. We isolate subsets of the entire sample, including low surface- 
gravity dwarfs, unusually blue L dwarfs, and photometric outliers in J — Kg color 
and investigate their kinematics. We find that the spectroscopically distinct 
class of unusually blue L dwarfs has kinematics clearly consistent with old age, 
implying that high surface- gravity and/or low metallicity may be relevant to their 
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spectral properties. The low surface-gravity dwarfs are kinematically younger 
than the overall population, and the kinematics of the red and blue ultracool 
dwarfs suggest ages that are younger and older than the full sample, respectively. 
We also present a reduced proper motion diagram at 2MASS Kg for the entire 
population and find that a limit of Hk^ > 18 excludes M dwarfs from the L and 
T dwarf population regardless of near-infrared color, potentially enabling the 
identification of the coldest brown dwarfs in the absence of color information. 

Subject headings: Astrometry- stars: low-mass- brown dwarfs- stars: fundamen- 
tal parameters 

1. INTRODUCTION 

Kinematic analyses of stars have played a fundamental role in shaping our picture of 
the Galaxy and its evolution. Prom early investigations (e.g. Schwarzschild 1908 , Lindblad 
1925, and Gort 1927) where the large scale structure of the Galactic disk was first explored, 
through more recent investigations (e.g. Gilmore & Reid 1983, Gilmore et al. 1989, Dehnen 
& Binney 1998, Famaey et al. 2005) where the structure of the Galaxy was refined to include 
a thick disk and prominent features such as streams, moving groups, and superclusters, 
kinematics have played a vital role in understanding the Galactic origin, evolution, and 
structure. Combining kinematics with spectral features, several groups have mapped out 
ages and metallicities for nearby F,G,K, and M stars (e.g. Nordstrom et al. 2004). The ages 
of these stars have become an important constraint on the Galactic star formation history 
and their kinematics have become a vital probe for investigating membership in the young 
thin disk, intermediate aged thick disk or older halo portion of the Galaxy. 

One population that has yet to have its kinematics exploited is the recently discovered 
population of very low-mass ultracool dwarfs (UCDs). These objects, which include those 
that do not support stable hydrogen fusion (Kumar 1962; Hayashi & Nakano 1963), oc- 
cupy the late- type M through T dwarf spectral classifications (e.g., Kirkpatrick 2005, and 
references therein). UCDs emit the majority of their light in the infrared and thus were 
only discovered in large numbers with the advent of wide-field near-infrared imaging sur- 
veys such as the Two Micron All Sky Survey (hereafter 2MASS; Skrutskie et al. 2006), the 
Deep Infrared Survey of the Southern Sky (hereafter DENIS; Epchtein et al. 1997) and the 
Sloan Digital Sky Survey (hereafter SDSS; York et al. 2000). Their very recent discovery has 
largely precluded astrometric measurements which require several-year baselines to produce 
useful measurements. Therefore, while UCDs appear to be comparable in number to stars 
(e.g., Reid et al. 1999), their role in the structure of the Galaxy has yet to be explored. 
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In addition, the thermal evolution of brown dwarfs (the lowest temperature ultracool 
dwarfs) implies that there is no direct correlation between spectral type and mass, leading 
to a mass/age degeneracy which makes it difficult to study the mass function and formation 
history of these objects. While some benchmark sources (e.g. cluster members, physical 
companions to bright stars) have independent age determinations, and spectroscopic analy- 
ses are beginning to enable individual mass and age constraints (e.g. Burgasser et al. 2006a; 
Saumon et al. 2007, Mohanty et al. 2004), the majority of brown dwarfs are not sufficiently 
characterized to break this degeneracy. Kinematics can be used as an alternate estimator 
for the age of the brown dwarf population. 

Moreover, kinematics can also be used to characterize subsets of UCDs. With hundreds 
of UCDs now known^, groupings of peculiar objects - sources whose photometric or spec- 
troscopic properties differ consistently from the majority of the population - are becoming 
distinguishable. Currently defined subgroups of late-type M, L and T dwarfs include 1) low 
surface- gravity, very low-mass objects (e.g. McCovern et al. 2004, Kirkpatrick et al. 2006, 
AUers et al. 2007, Cruz et al. 2007), 2) old, metal-poor ultracool subdwarfs (e.g. Burgasser 
et al. 2003a, Lepine et al. 2003, Gizis & Harvin 2006, Burgasser et al. 2007b), 3) unusually 
blue L dwarfs (hereafter UBLs; e.g. Cruz et al. 2003, Cruz et al. 2007; Knapp et al. 2004, 
Chiu et al. 2006), and 4) unusually red and possibly dusty L dwarfs (e.g. Looper et al. 2008; 
McLean et al. 2003). While observational peculiarities can overlap between these groups 
(e.g. both young and dusty L dwarfs can be unusually red), they appear to encompass ob- 
jects with distinct physical traits (e.g., mass, age, composition, and cloud properties) so they 
are important for drawing a connection between observational characteristics and intrinsic 
physical properties. Kinematics can be used to investigate the underlying physical causes 
for the peculiarities of these groups. 

In the past decade a number of groups have conducted astrometric surveys of UCDs, 
including subsets of low mass objects (e.g. Vrba et al. 2004, Dahn et al. 2002, Gizis et al. 
2000, Tinney et al. 2003, Schmidt et al. 2007, Jameson et al. 2007, Osorio et al. 2007, and 
West et al. 2008, 2006). We have initiated the Brown Dwarf Kinematics Project (BDKP) 
which aims to measure the positions and 3D velocities of all known L and T dwarfs within 20 
pc of the Sun and selected sources of scientific interest at larger distances (e.g. low surface- 
gravity dwarfs, subdwarfs). In this article we add 332 new proper motion measurements and 
combine all published proper motion measurements and distance estimates into a uniform 
sample to examine the ultracool dwarf population as a whole. Section 2 of this paper outlines 
the observed sample and describes how proper motion measurements were made. Section 3 
discusses the expanded sample and how distances and Vfan measurements were calculated. 



^ An up-to-date list of known L and T dwarfs is maintained by C. Gelino, D. Kirkpatrick and A. Burgasser 
at http://www.dwarfarchives.org 
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Section 4 examines the full astrometric sample and subsets. Section 5 reviews the high 
tangential velocity objects in detail. Finally, section 6 applies an age-velocity relation and 
examines resultant ages of the full sample and red/blue outliers. 

2. OBSERVATIONS AND PROPER MOTION MEASUREMENTS 

2.1. Sample Selection 

Our goal is to re-image all known late-type M, L, and T dwarfs to obtain accurate 
uniformly measured proper motions for the entire ultracool dwarf population. In our sample 
we focused on the lowest temperature L and T dwarfs that were lacking proper motion 
measurements or whose proper motion uncertainty was larger than 40 mas/yr. We gave high 
priority to any dwarf that was identified as a low surface-gravity object in the literature. Our 
sample was created from 634 L and T dwarfs listed on the Dwarf Archives website as well as 
456 M7-M9.5 dwarfs gathered from the literature (primarily from Cruz et al. 2003, 2007). 
The sample stayed current with the Dwarf Archives website through April 2008. Figure 1 
shows the histogram of spectral type distributions for the entire sample. The late-type M 
and early-type L dwarfs clearly dominate the ultracool dwarf population. Plotted in that 
figure is the current distribution of objects with proper motion values and the distribution of 
objects for which we report new proper motions. To date we have re-imaged 427 objects. As 
of June 2008 and including all of the measurements reported in this article, 570 of the 634 
known L and T dwarfs and 277 of the 456 late-type M dwarfs in our sample have measured 
proper motions. 

2.2. Data Acquisition and Reduction 

Images for our program were obtained using three different instruments and telescopes 
in the northern and southern hemispheres. Table 1 lists the instrument properties. For 
the northern targets the 1.3m telescope at MDM with the TIFKAM IR imager in J band 
was used. For the southern targets the 0.9m and 1.5m telescopes at CTIO with the CFIM 
optical imager in I band and the CPAPIR wide field-IR imager in J band (respectively) 
were used. The CTIO data were acquired through queue observing on 11 nights in March, 
September and December 2007, and standard user observing on 9 nights in January 2008. 
The MDM targets were imaged on five nights in November 2007 and 7 nights in April 2008. 
Objects were observed as close to the meridian as possible up to an airmass of 1.80, and with 
seeing no greater than 2.5" FWHM. Exposure times varied depending on the target and the 
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instrument. For CPAPIR the exposure times ranged over 15— 40s with 4 co-adds per image 
and a five-point dither pattern. At MDM the exposure times ranged over 30— 120s with up 
to 6 co-adds per image and a three to five-point dither pattern. For the 0.9m observations 
the exposure times ranged over 180— 1800s per image with no co-adds and a three-point 
dither pattern. The dither offset between positions with each instrument was 10". 

All data were processed in a similar manner using standard IRAF and IDL routines. 
Domcflats were constructed in the J or I band. CPAPIR and CFIM domcflats were cre- 
ated from 10 images illuminated by dome lamps, and TIFKAM domcflats were created by 
subtracting the median of 10 images taken with all dome lights off from the median of 10 
images taken with the dome lights on. A dark image constructed from 25 images taken with 
the shutter closed was used to map the bad pixels on the detector. Domefiats were then 
dark subtracted and normalized. Sky frames were created for each instrument by median- 
combining all of the science data taken on a given night. Science frames were first fiat-fielded, 
then sky-subtracted. Individual frames were shifted and stacked to form the final combined 
images. 



2.3. Calculating Proper Motions 

The reduced science frames were astrometrically calibrated using the 2MASS Point 
Source catalogue. 2MASS astrometry is tied to the TYCHO 2 positions and the reported 
astrometric accuracy varies from source to source. In general the positions of 2MASS sources 
in the magnitude range 9 < < 14 are repeatable to 40-50 mas in both RA and DEC. 

Initial astrometry was fit by inputting a 2x2 transformation matrix containing astrom- 
etry parameters which were first calculated from one image in which two stars with known 
2MASS RA and DEC values and second epoch X,Y pixel positions were known. The refer- 
ence RA, DEC and pixel values were first set to the pointing RA and DEC values and the 
center of the chip respectively. 

RA and DEC values for all stars in the field were then imported from the 2MASS 
point source catalogue and converted to X,Y pixel positions using the initial astrometric 
parameters. We worked in X,Y positions of the second epoch image so we could overplot 
point source positions on an image and visually check that we converged upon a best fit 
solution. We detected point sources on the second epoch image with a centroiding routine 
which used a detection threshold of 5 sigma above the background. We matched the 2MASS 
X,Y positions to the second epoch positions by cross-correlating the two lists. We refined 
the astrometric solution by a basic six parameter, least-squares, linear transformation where 
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we took the positions from the 2MASS image {Xi,Yi) and the positions from the second 
epoch image {Xq,Yq) and solved for the new X,Y pixel positions of the second epoch image 
in the 2MASS frame. Due to the large field of view, we checked for higher order terms in 
the CPAPIR images and found no significant terms. The following equations were used: 

X^X2o + A(X, - Xo) + B(Y, - Yo) (1) 

Y = y2o + C{Xi - Xo) + D{Yi - Fq) (2) 

where X2o and were set to the center of the field, A,B,C,D solve for the rotation and 
plate scale in the two coordinates. 

The sample of stars used to compute the astrometric solution for each image were 
selected according to the following criteria: 

1. Only stars in the 2MASS J magnitude range 12 < J < 15 were used, as objects in this 
intermediate magnitude range transformed with the smallest residuals from epoch to 
epoch. 

2. The solution reference stars were required to transform with total absolute residuals 
against 2MASS of < 0.2 pixels. From testing with images taken consecutively using 
each instrument, the best astrometric solution always generated between 0.1 and 0.2 
pixel average residuals. Therefore the stars used to calculate the solution were required 
to fall in or below that range. 

As the solution was iterated, the residuals were examined at each step, and stars that 
did not fit the above criteria were removed. For CPAPIR, the process converged on a solution 
that had between 100 and 200 reference stars with average residuals < 0.15 pixels. TIFKAM 
and CFIM have smaller fields of view (~6 arcmin and 5 arcmin respectively as opposed 
to 35 arcmin for CPAPIR) so there were far fewer stars to work with. For these imagers the 
process converged on a solution that had between 15 and 60 reference stars. The astrometric 
solution was required to converge with no less than 15 reference stars and when this criterion 
could not be met, the other two criteria above were relaxed. As a result TIFKAM and CFIM 
had shghtly larger residuals on the astrometric solution (average residuals < 0.25 pixels). 

Once an astrometric solution was calculated, final second epoch positions were computed 
using a Gaussian fit for each 2MASS X,Y position on an image. For the science target, a 
visual check was employed to ensure that it had been detected and X,Y positions were 
manually input for the Gaussian fit. Final X,Y positions were then converted back into RA 
and DEC values using the best astrometric solution and the proper motion was calculated 
using the positional offset and time difference between the second epoch image and 2MASS. 
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The residuals of the astrometric solution were converted into proper motion uncertainties 
by multiplying by the plate scale of the instrument and then dividing by the epoch difference. 
The baselines ranged from 6-10 years and our astrometric uncertainties range from 5 to 50 
mas/yr. Positional uncertainties for each source were also calculated by comparing the 
residuals of transforming the X,Y position for our target over consecutive dithered images. 
These uncertainties are dominated by counting statistics with the high S/N sources having 
negligible positional uncertainties compared to the uncertainties in the astrometric solution. 
We added the positional and astrometric solution uncertainties in quadrature to determine 
the total proper motion uncertainty. Figure 2 shows the distribution of proper motion 
uncertainties and baselines for all new proper motion measurements reported in this paper. 
The median uncertainty was 18 mas/yr. 

Of the 427 proper motion measurements we report in this paper, 332 are presented here 
for the first time. Twelve objects were purposely remeasured with multiple instruments as a 
double check on the accuracy of the astrometric solution, and 42 objects were remeasured to 
refine the proper motion uncertainties. Thirty-two measurements were published in Jameson 
et al. (2007)-hereafter J07-, and 11 in Caballero (2007) while our observations were underway. 
The proper motion measurements presented in this paper agree to better than 2a in 84 of 
the 97 cases of objects with prior measurements. Table 4 lists those cases where the proper 
motions are discrepant by more than 2a with a pubhshed value. For nine of the objects, 
there is a third (fourth or fifth) measurement by an independent group with which we are 
in good agreement. We are discrepant with six objects reported in Deacon et al. (2005) 
but we note that there are no position angle uncertainties reported for these objects in that 
catalogue therefore we cannot fully assess the accuracy of the proper motion components. 
The difference in proper motion for 2MASSW J1555157-095605 is quite large (> 1" /yi 
difference) but there are two other measurements for this object with which we are in close 
agreement. We have examined all of the discrepant proper motion images carefully and see 
no artifacts that could have skewed our measurements. Figure 3 compares the proper motion 
component measurements from this paper with those from the literature for objects with 
< 0.5 "/yr and jierr < 0.1 "/y^- With~90% agreement with published results, this indicates 
that the 332 new measurements are robust. Table 2. contains all new measurements reported 
in this article, and Table 3. contains the astrometric measurements for the full sample. 
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3. DISTANCES, TANGENTIAL VELOCITIES, AND REDUCED PROPER 

MOTION 

3.1. Expanded Sample 

We extended our observational sample to include published late-type M, L, and T dwarfs 
with proper motion measurements yielding a full combined sample containing 841 objects. 
Thirty-three percent of ultracool dwarfs in the full sample have multiple proper motion 
measurements. In these cases we chose the measurement with the smallest uncertainty for 
our kinematic analysis, typically objects from high precision astrometric surveys such as Vrba 
et al 2004 or Dahn et al 2002. If there was a value discrepant by more than 2a amongst 
multiple measurements (> 2) for an object then regardless of uncertainty we defaulted to the 
numbers that were in agreement and chose the one with the smaller uncertainty. Otherwise, 
if there was a discrepancy and only two measurements, we quoted the one that had the 
smaller uncertainty and made note of it during the analysis. 



3.2. Distances and Tangential Velocities 

True space velocities arc a more fundamental measure of an object's kinematics than 
apparent angular motions, so proper motions for the complete sample were converted to 
tangential velocities using astrometric or spectrophotometric distances. As of January 2008, 
only 79 of the 634 L and T dwarfs and 64 of the 456 late-type M dwarfs in our sample had 
published parallax measurements. Therefore to include the other 83% of L and T dwarfs and 
87% of late- type M dwarfs in a population analysis, published absolute magnitude/spectral 
type relations were used for calibrating distances. Dahn et al. (2002) and Vrba et al. (2004) 
both showed that Mj is well correlated with spectral type for late-type M, L, and T dwarfs 
(see also West et al. 2005, and Covey et al. 2007). Since the initial relations were published 
several investigators have revised the absolute magnitude/spectral type relation after in- 
cluding new measurements and removing resolved binaries. In this paper, the distances for 
the M7-L4.5 dwarfs were calculated using the absolute 2MASS J magnitude/spectral type 
relation in Cruz et al. (2003) and the distances for the L5 -T8 dwarfs were calculated using 
the absolute MKO K magnitude/spectral type relation in Burgasser (2007)^. Both optical 



^The coefficients of this polynomial relation reported in Burgasser (2007) did not list sufficient significant 
digits, yielding a slightly different numerical relation than that used in the paper's analysis. The coefficients 
as defined should be {q} = [10.4458, 0.232154, 0.0512942, -0.0402365, 0.0141398, -0.00227108, 0.000180674, 
-6.985016-06, 1.051196-07], where Mk = ELoCiSpT' and SpT(TO) = 10, SpT(T5) = 15, etc. 
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and near-IR spectral types are reported for ultracool dwarfs. For late-type M through the 
L dwarfs, we use the optical spectral type in the distance relation when available but use 
near-IR spectral types when no optical spectral types are reported. We use the near-IR 
spectral type in the distance relation for all T dwarfs. The Cruz et al. (2003) relation was 
derived for the 2MASS magnitude system, while the Burgasser (2007) relation was derived 
using the MKO system. In reporting distances we maintain the magnitude system for which 
the relation was calculated, converting a 2MASS magnitude to an MKO magnitude or vice 
versa using the relation in Stephens & Leggett (2004) when necessary. The most recent 
precision photometry for many L and T dwarfs (e.g. Knapp et al 2004; Chiu et al. 2006, 
2008) are reported on the MKO system; yet the majority of objects explored in this paper 
have measured 2MASS magnitudes. We convert MKO filter measurements to the 2MASS 
system when available using the conversion relations of Stephens & Leggett (2004) so all of 
the ultracool dwarf photometry in Table 3 is reported on the 2MASS system. 

The uncertainty in the derived distance is dominated by the uncertainty in the spectral 
type (the photometric uncertainties are typically between 0.02-0.1 mag whereas the spectral 
type uncertainties are typically 0.5-1.0). This leads to a systematic over- or underestima- 
tion of distance of up to 30%. Therefore the kinematic results presented in this paper are 
largely sensitive to the reliability of the spectrophotometric distances used to calculate Vtan- 
Furthermore, unresolved multiplicity leads to an underestimation of distance. Recent work 
has shown that roughly 20% of ultracool dwarfs are likely to be binary (Allen 2007, Reid 
et al. 2008), and this fraction may be even higher across the L dwarf/T dwarf transition 
(Burgasser et al. 2006b). Seven percent (56) of the dwarfs analyzed in this paper are known 
to be close binaries and of these, most appear to be near equal-mass/equal brightness (e.g. 
Bouy et al. 2003; Burgasser et al. 2006b). For these objects we use the distances quoted in 
the binary discovery papers where the contribution of flux from the secondary was included 
in the distance estimate. Any remaining tight binaries probably constitute no more than 
10-20% of the sample and the contamination of their inclusion in the kinematic analysis is 
relatively small. 



3.3. Reduced Proper Motion Diagram 

A reduced proper motion diagram is a useful tool for distinguishing between kinematically- 
distinct stellar and substellar populations. This parameter was used extensively in early high 
proper motion catalogues to explore Galactic structure (Luyten 1973). Proper motion is 
used as a proxy for distance measurements following the expectation that objects with large 
proper motions will be nearest to the Sun. The definition is analogous to that of absolute 
magnitude: 



if = m + 5.0 + d.Ologioiii) 



(3) 



or 

H^M + 5.0logio{Vtan) - 3.38 (4) 

where m and M are the apparent and absolute magnitudes (respectively), Vtan is measured 
in km/s and /i is measured in "/yr. 

We can use reduced proper motion to search for the lowest temperature objects. In 
Figure 4, we show the reduced proper motion at Kg for our astrometric sample. We find 
that below an H^^ of 18 there are only L and T dwarfs regardless of near-IR color. Since the 
discovery of the first brown dwarfs, near-IR color selection has been the primary technique 
for identifying strong candidates. But because M dwarfs dominate photometric surveys (they 
are bright, nearby and found in large numbers), near-IR color cut-offs were administered to 
maximize the L and T dwarfs found in searches. These cut-offs have caused a gap in the 
near-IR color distribution of the brown dwarf population, particularly around J — i^T^ of 1 
where early-type T dwarfs and metal weak L dwarfs are eliminated along with M dwarfs. 
A reduced proper motion diagram with the cut-off limit cited above allows a search which 
eliminates the abundant M dwarfs and probes the entire range oi J — Kg colors for the 
ultracool dwarf population. 

Note that, while our cut-off limits are good guidelines for segregating the coolest tem- 
perature dwarfs within the ultracool dwarf population, there is likely to be contamination 
in selected regions of the sky from relatively rare ultracool subdwarfs and cool white dwarfs, 
which are nonetheless of scientific interest. 



4. ANALYSIS 

4.1. Kinematic Characteristics of the Ultracool Dwarf Population 

The ultracool dwarfs analyzed in this paper have a range of proper motion values from 
0.01— 4.7"/yr and a range of proper motion uncertainties from 0.0002 — 0.3 " /jr. While one 
of our goals is to refine proper motion measurements of ultracool dwarfs to have uncertainties 
< 40 mas/yr, there arc still 86 or 10% that have larger errors. Since the uncertainty in Vtan 
is generally dominated by the uncertainty in distance (see subsection 3.2 above) we make 
no restrictions on the accuracy of the proper motion measurements used in the kinematic 
analysis. The median la detection limit for proper motion measurements in this paper was 
18 mas/yr (see Figure 2). We use this number as a proxy for the L and T dwarfs (where 
we are looking at all known field objects as opposed to the late-type M dwarfs where we are 
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looking at only a subset) to determine the percentage of objects with appreciable motion. 
We find that 32 move slower than our 2(T detection limit and ten of those are at or below 
our la limit. This indicates that according to our astrometric standard, less than 6% of L 
and T dwarfs have no appreciable motion. Conversely, 32 objects (or 6% of the population) 
move faster than 1.0"/yr making them some of the fastest proper motion objects known. As 
late-type dwarfs are intrinsically quite faint and have only been detected at nearby distances 
(generally < 60 pc), the high proper motion values measured are not surprising. Using the 
median proper motion values listed in Table 5 as a proxy , we can also conclude that at 
least half or more of the brown dwarf population would be easily detectable on a near-IR 
equivalent of Luyten's 2-tenth catalog (Luyten 1979) where the limiting proper motion was 
~ 0.15 "/yr. 

Table 5 lists the average proper motion values and photometric data for the entire 
population binned by spectral type. There is a trend within these data for larger proper 
motion values with increasing spectral type. This is clearest within the L0-L9 population 
where the sample is largest. We further bin this group into thirds to compare a statistically 
significant sample. We examine the L0-L2, L3-L5, and L6-L9 populations and find the 
median proper motion values increase as (0.174, 0.223, 0.289) "/jr. This trend most likely 
reflects the fact that earlier type sources are detected to further distances. Indeed when 
we examine the median distance values for these same groupings we find values of (31, 27, 
20) pc respectively. 



4.2. Kinematics of Full and 20pc Samples 

We have conducted our kinematic analysis on two samples: the full astrometric sample 
and the 20 pc sample. Figure 5 shows the distance distribution for all ultracool dwarfs 
regardless of proper motion measurements to demonstrate the pertinence of the 20 pc sample. 
In this figure, both the late-type M and L dwarfs diverge from an N oc d^ density distribution 
around 20 pc. The T dwarfs diverge closer to 15 pc. Within the hterature (e.g Cruz et al. 
2003) complete samples to 20 pc have been reported through mid-type L dwarfs so we use 
this distance in order to establish a volume-limited kinematic sample. Wc also examine the 
two samples with and without objects with Vtan > 100 km s~^ in order to remove extreme 
outliers that may comprise a different population. 

Tables 6-7 contain the mean kinematic properties for the 20 pc sample and the full 
astrometric sample respectively. Figure 6 shows Vtan vs. spectral type for both samples. 
As demonstrated in Figure 6, we find no difference between the two samples, with median 
Vtan values of 26 km s~^ and 29 km s~^ and atan values of 23 km s~^ and 25 km s~^ for 
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M7-T9 within the full sample and the 20 pc sample respectively. Within spectral class bins, 
namely the M7-M9, L0-L9, or T0-T9 groupings, we find no significant kinematic difi^erences. 
This indicates that we are sampling a single kinematic population regardless of distance and 
spectral type. 

Figure 7 shows the distribution of tangential velocities. There are 14 objects with Vtan 
> 100 km that fall at the far end of the distribution. Exclusion of these high velocity 
dwarfs naturally reduces the median Vtan and atan values. The most significant difference in 
their exclusion occurs within the L0-L9 group as 10 of the 14 objects belong to that spectral 
class. We explore the importance of this subset of the ultracool dwarf population in section 5. 

In order to put our kinematic measurements in the context of the Galaxy we compare 
with Galactic U,V,W dispersions. Proper motion, distance, and radial velocity are all re- 
quired to compute these space velocities. Therefore, a direct Galactic U,V,W comparison 
with the ultracool dwarf population is not possible because radial velocity measurements for 
ultracool dwarfs are sparse, with only 48 of the L and T dwarfs to date having been reported 
in the literature (e.g. Mohanty & Basri 2003, Osorio et al. 2007, Bailer- Jones 2004). This 
is a similar problem to that for precise brown dwarf parallax measurements, but there is no 
relationship for estimating radial velocities as there is for estimating distances. However, we 
can divide our sample into three groups along Galactiocentric coordinate axes (toward poles, 
in direction of Galactic rotation and radially to/from Galactic center) in order to minimize 
the importance of radial velocity in 2 out of the 3 space velocity components. Wc create 
cones of (all inclusive), 30, and 60 degrees around the galactic X,Y, and Z axes. Inside of 
each cone we set either the U,V, or W velocity to zero if the cone surrounds the galactic X,Y 
or Z axis respectively. In this way we can set the radial velocity of each source to zero with 
minimum impact on the component velocities of the entire sample and gather U,V,W infor- 
mation for the known ultracool dwarf population. We emphasize that this analysis is crude 
as the distribution of ultracool dwarfs is not isotropic (the Galactic plane has largely not been 
explored), and while the cones help to minimize the importance of radial velocity unless an 
object is directly on the X,Y, or Z axis, the radial velocity component will contribute to the 
overall velocities. Therefore, the spread of U,V,W velocities will be biased toward a tighter 
dispersion than the true values. In order to calculate total velocities (Vtot) for objects, which 
requires U,V, and W velocities we choose a cone of 30 degrees which provides a statistically 
significant sample. We create the cone around the X,Y, or Z axis and assume that within 
that cone either the V,W or U,Z or U,V components respectively are correct. To obtain the 
third component we assume it to be the average of the two calculated ones. In this way we 
can gather Vtot information which will be used for age calculation purposes in Section 6. 

Figure 8 shows our resultant U,V,W distributions where we measure {au, (Jy, (7vf) = (28, 
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22, 17) km s~^. We compare these dispersions with the kinematic signatures of the three 
Galactic populations, namely the thin disk, the thick disk, and the halo. The overwhelming 
majority of stars in the solar neighborhood are members of the Galactic disk and these are 
primarily young thin disk objects as opposed to older thick disk objects. The halo population 
of the Galaxy encompasses the oldest population of stars in the Galaxy but these objects 
are relatively sparse in the vicinity of the Sun. Membership in any Galactic population has 
implications for the age and metallicity of the object and kinematics play a large part in 
defining the various populations. Soubiran et al. (2003) find (cr;/, cry, (T|^)= (39 ± 2, 20 ± 2, 
20 ± 1) km s-^ for the thin disk and {uu, oy, (7vk)=(63 ± 6, 39 ± 4, 39 ± 4) km s^^ for the 
thick disk, and Chiba & Beers (2000) find (cTf/, ay, (7w^)=(141 ± 11, 106 ± 9, 94 ± 8) km s"^ 
for the halo portion of the Galaxy. Our U,V,W dispersions are consistent (albeit narrower 
in U ) with the Galactic thin disk. 

Osorio et al. (2007) - hereafter Os07- examined 21 L and T dwarfs and found [ou, ay, (7vt^)=(30.2, 
16.5, 15.8) km s~^. Their velocity dispersions are tighter than what is expected from the 
Galactic thin disk population. Our calculated dispersions are tighter at U than the Os07 
result (which is expected due to the stated bias) but broader in V and W. In section 6 we 
discuss the implications on age of the differences calculated from our astrometric sample. 

4.3. Red and Blue Photometric Outliers 

As discussed in Kirkpatrick et al. (2005) the large number of late-type M, L, and T 
dwarfs discovered to date has revealed a broad diversity of colors and spectral characteris- 
tics, including specific subgroups of peculiar sources that are likely related by their common 
physical properties. As a very basic metric, near-IR colors provide one means of distinguish- 
ing between "normal" and "unusual" objects. To investigate our sample for kinematically 
distinct photometric outliers, we first defined the average color ((J — Ks)avg) ^ well as stan- 
dard deviation (aj^Ks) ^ ^ function of spectral type using all known ultracool dwarfs (i.e., 
both with and without proper motion measurements). Defining the {J — Ks)avg for spectral 
bins has been done in previous ultracool dwarf studies such as Kirkpatrick et al. (2000), 
Vrba et al. (2004), and West et al. (2008) but wc have included all ultracool dwarfs in the 
dwarfarchives complilation and the updated photometry reported in Chiu et al. (2006, 2008) 
and Knapp et al. (2004) which we have converted from the MKO system to the 2MASS 
system. Objects were eliminated from the photometric sample if they fit any of the following 
criteria: 

1. Uncertainty in ^ or Kg > 0.2 magnitude; 
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2. Known subdwarf; 

3. Known binaries unresolved by wide-field imaging surveys (i.e. separations ^ 1 " e.g. 
Martin et al. 1999; Bouy et al. 2003; Burgasser et al. 2006b; Close et al. 2003; Liu 
et al. 2006; Reid et al. 2006); and 

4. Member of a star forming region (such as Orion) or open cluster (such as the Pleiades) 
indicating an age < 100 Myr (e.g. AUers et al. 2007; Zapatero Osorio et al. 2002) 

We then designated objects as photometric outliers if they satisfied the following criterion : 
Aj-i^, =1 (J - X,) - (J - K,)a,g \> max(2aj_x., 0.4) (5) 

In other words, if an object's J — Kg color was more than twice the standard deviation of 
the color range for that spectral bin than we flagged it as a red or blue photometric outlier. 
If twice the standard deviation was larger than 0.4 magnitudes then it was automatically 
reset to 0.4. We chose 0.4 as the maximum upper limit for 2aj-Ks this is the Aj^k^ for 
the entire ultracool dwarf population. 

There are relatively few objects in each spectral bin beyond L9. For spectral type (SpT) 
< L9 there is a mean of 45 objects used per bin whereas for SpT> L9 there is a mean of only 
7 objects. So photometric outliers arc difficult to define for the lower temperature classes 
and may contaminate the analysis. We grouped T7-T8 dwarfs to improve the statistics used 
to calculate average values. The kinematic results for this subset of the ultracool dwarf 
population are reported with and without the T dwarfs in Table 8. Figure 9 shows the 
resulting J — Kg color distribution and highlights the photometric outliers. Tables 9-10 list 
the details of the red and blue photometric outliers respectively. Table 5 lists the resultant 
mean photometric values for each spectral type. 

Amongst the full sample, we find 16 blue photometric outliers and 29 red photometric 
outliers. Many of the objects have already been noted in the literature as having unusual 
colors, and several of these have anomalous spectra and have been analyzed in detail (e.g. 
Burgasser et al. 2007a, Knapp et al. 2004, Folkcs et al. 2007, Chiu et al. 2006). Table 8 
lists the mean kinematic properties for the blue and red subgroups of the ultracool dwarf 
population and Figure 10 isolates the outliers and plots their tangential velocity vs. spectral 
type. The blue outhers have a median Vtan value of 53 km and a atan of 47 km s~^ while 
the red outliers have a median Vtan value of 26 km s~^ and a atan of 16 km s~^ . Figure 
11 shows the tangential velocity vs. 3- Kg deviation for all objects in the sample with the 
dispersions of the red and blue outhers highlighted. There is a clear trend for Vtan values 
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to decrease from objects that are blue for their spectral type to those that are red. This is 
particularly significant at the extreme edges of this diagram. The dashed line in FigTirc 1 1 
marks the spread of Vtan values for the full sample and demonstrates the significant deviations 
for the color outliers. We explore the age differences from these measurements in section 6. 

4.4. Low Gravity Objects 

A number of ultracool dwarfs that exhibit low surface-gravity features have been re- 
ported in the literature within the past few years (e.g. Cruz et al. 2007, Luhman & Rieke 
1999; McGovern et al. 2004, Kirkpatrick et al. 2006; AUers et al. 2007). Low surface-gravity 
dwarfs are distinguished as such by the presence of weak alkali spectral features, enhanced 
metal oxide absorption, and reduced H2 absorption. They arc most likely young , with lower 
masses than older objects of the same spectral type. For ages < 100 Myr these objects may 
also have larger radii than older brown dwarfs and low mass stars with similar spectral types, 
as they are still contracting to their final radii (e.g.. Burrows et al. 1997). 

We examine the kinematics of 37 low surface-gravity dwarfs in this paper. Seven of these 
objects are flagged as red photometric outhers and were examined in the previous subsection. 
The overlap between these two subgroups is not surprising as the reduced H2 absorption in 
low surface-gravity dwarfs leads to a redder near-IR color. The median Vtan value for this 
subgroup is 18 km and the atan value is 15 km s~^ which is smaller than that of the red 
photometric outliers as a whole and therefore points to the same conclusion. The smaller 
median Vtan and tighter dispersion of the low surface-gravity dwarfs as compared to either 
the full or 20 pc sample indicates that they are kinematically distinct. 

4.5. Unusually Blue L Dwarfs 

A subgroup of unusually blue L dwarfs (UBLs) has been distinguished based on strong 
near-IR H2O, FeH and K I spectral features but otherwise normal optical spectra. Burgasscr 
et al. (2008a) -hereafter B08- identify ten objects that comprise this subgroup (see Table 6 
in B08). With the kinematics reported in this article we are able to analyze all ten. There are 
several physical mechanisms that can contribute to the spectral properties of UBLs. High 
surface-gravity, low metallicity, thin clouds or unaccounted multiplicity are amongst the 
physical mechanisms most often cited. B08 has demonstrated that while sub-solar metaUicity 
and high surface-gravity could be contributing factors in explaining the spectral deviations, 
thin, patchy or large-grained condensate clouds at the photosphere appears to be the primary 
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cause for the anomalous near-IR spectra (e.g. Ackerman & Marley 2001, Burrows et al. 2006). 

The median Vtan value for this subgroup is 99 km s~-^ with atan of 47 km s"-*^ and this 
subgroup consists of dwarfs with the largest Vtan values measured in this kinematic study. 
These kinematic results strengthen the case that the UBLs represent an older population 
and that the blue near-IR colors and spectroscopic properties of these objects are influenced 

by large surface-gravity and/or slightly subsolar metallicities. Both of these effects may be 
underlying explanations for the thin clouds seen in blue L dwarf photospheres. Subsolar 
metallicity reduces the elemental reservoir for condensate grains while high surface-gravity 
may enhance gravitational settling of clouds. In effect, the clouds of L dwarfs may be tracers 
of their age and/or metallicity. 

Eight of the ten UBLs examined in this subsection are also flagged as blue photometric 
outliers and examined in detail above. The overlap between these two subgroups is not 
surprising as many of the UBLs were initially identified by their blue near-IR color (e.g. 
Cruz et al. 2007, Knapp et al. 2004). There are 8 other blue photometric outliers, one of 
which has a Vtan value exceeding 100 km s~^. We plan on obtaining near-IR spectra for these 
outliers to investigate the possibility that they exhibit similar near-IR spectral features to 
the UBLs. 

While the UBLs are the most kinematically distinct subgroup analyzed in this paper, 
their kinematics do not match those of the ultracool subdwarfs. The subdwarfs were excluded 
from the kinematic analysis in this paper because they are confirmed members of a separate 
population. The median Vtan value for this subgroup is 196 km s~^ with atan of 91 km s~^. 
The UBLs move at half this speed indicating there is a further distinction between UBLs 
and the metal-poor halo population of ultracool dwarfs. 

5. High Velocity Dwarfs 

Table 11 summarizes the properties of the 14 high velocity dwarfs whose Vtan measure- 
ments exceed 100 km s~^. A number of these have been discussed in the literature, having 
been singled out in their corresponding discovery papers as potential members of the thick 
disk or halo population. One high velocity dwarf is presented here for the first time. SDSS 
J093109. 56+032732. 5 is an L7.5 dwarf and is classified as both a UBL and a blue photometric 
outlier. We calculate Vtan for this object to be 108 ± 23 km s~^. 

Among the high velocity dwarfs, 11 have colors that are blue and 3 have colors that are 
normal for their spectral type. Three objects belong to the UBL subgroup. Three of the 
objects are late-type M dwarfs (2MASS J18261131+3014201, 2MASS J03341218-4953322, 
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and 2MASS J132352+301433), one is a late T7.5 dwarf (2MASSJ 11145133-2618235) and 
the rest are early to mid-type L dwarfs. Four of the objects are flagged as blue photometric 
outliers. We explore the possibility that these objects are thick disk or halo objects in detail 
in a forthcoming paper. 



6. ON THE AGES OF THE ULTRACOOL DWARF POPULATIONS 

6.1. Kinematics and Ages 

A comparison of the velocity dispersion for nearby stellar populations can be an indicator 
of age. While individual Vtan measurements cannot provide individual age determinations 
due to scatter and projection effects, the random motions of a population of disk stars are 
known to increase with age. This effect is known as the disk age-velocity relation (AYR) 
and is simulated by fitting well-constrained data against the following analytic form: 

a{t)^ao{l + -r (6) 

T 

where a{t) is the total velocity dispersion as a function of time, o"o is the initial velocity 
dispersion at t=0, r is a constant with unit of time, and a is the heating index (Wielen 1977). 
a{t) is defined for U,V,W space velocities but we can estimate the total velocity dispersion 
using our measured tangential velocities assuming the dispersions are spread equally between 
all three velocity components, such that: 

a{t) = (3/2)(V2)^^^^ 

Hanninen & Flynn (2002) calculate a from seven distinct data sets (both pre and post- 
Hipparcos) and find a ranges from 0.3 to 0.6. This is a large range of values and the authors 
are reluctant to assign a higher likelihood to any given value as each have nearly equal 
uncertainties. One possible explanation for the spread of values is that a should be mass 
dependent ^. If so, this would make a large difference in the age calculations for the low 
mass ultracool dwarf population. While the age- velocity relation in the nearby disk remains 
only roughly determined, there is strong observational evidence for a relation so we proceed 
with caution in examining the broad age possibilities implied by the AVR for the ultracool 
dwarf population. 



^Indeed Iwanowska 1980 proposes the introduction of a mass term to account for the importance of the 
exchange of energy between stars in the Galactic disk 



- 18 - 



Recent findings have suggested that late-type M stars in the solar neighborhood are 
younger on average than earlier type stars (Hawkins & Bessell 1988; Kirkpatrick et al. 1994; 
Reid et al. 1994). Several investigators have combined kinematics with the Wielen (1977) 
relationship (which uses a value of 1/3 for a) to estimate age ranges for the ultracool dwarf 
population and concluded that it is kinematically younger than nearby stellar populations 
(e.g. Dahn et al. 2002, Schmidt et al. 2007, Gizis et al. 2000, Osorio et al. 2007). We 
conducted a direct Vtan comparison with nearby stellar populations to draw conclusions 
about the kinematic distinguishability of our ultracool dwarf sample. We compared the 
kinematics of a 20 pc sample of F,G,K, and early M stars from Soubiran et al. (2003), 
Kharchenko et al. (2004) and Nordstrom et al. (2004) using a limiting proper motion of 25 
mas/yr to our 20 pc sample and examined the resultant median Vtan, (^tan, and Vtot , o^tot 
values (where Vtot comes from the U,V,W velocities). Figure 12 shows our resultant velocity 
dispersions for nearby stellar populations along with the dispersions of our 20 pc sample. We 
show both the dispersions calculated using tangential velocities and those calculated using 
U,V,W velocities. As expected the dispersions are tighter for the UCDs when U,V,W values 
are used since we have attempted to minimize the importance of radial velocity. This effect is 
also reflected in the younger ages estimated from these dispersions. The tangential velocity 
dispersions are in good agreement between the UCDs and nearby stellar populations (see 
also West et al. 2008,Bochanski et al. 2007, and Covey et al. (2008)). Table 12 contains the 
calculated kinematic measurements and Wielen ages using both Utan and atot- With Wielen 
ages of 3-8 Gyr calculated from atan-, we conclude that our 20 pc sample is kinematically 
indistinct from other nearby stellar populations and hence is not kinematically younger. The 
ages calculated by the AVR for the 20 pc sample are in good agreement with those predicted 
in population synthesis models where the mean ages for the ultracool dwarf population range 
from 3-6 Gyr (Burgasser 2004b; Allen et al. 2005). 

We do find younger ages for the ultracool dwarf population when the high velocity 
dwarfs are excluded. As stated in section 4, the median Vtan and atan values are naturally 
reduced when the high velocity dwarfs are excluded and consequently the ages are also 
reduced. Kinematic analyses of the past have regarded these objects as a separate older 
population and omitted them from an age calculation (e.g. Schmidt et al. 2007). Table 6 
presents the ages with and without the high velocity dwarfs for the 20 pc sample. When 
the high velocity dwarfs are excluded the age ranges are reduced from 3-8 Gyr to 2-4 Gyr, 
which is still consistent with population synthesis models. 

The Os07 study estimated mean ages of ~ 1 Gyr for the L and T dwarf population . Even 
with the exclusion of the kinematic outliers, the ages calculated in our full and 20 pc samples 
do not match this very young age. Os07 combined proper motions, precise parallaxes, and 
radial velocities to study the 3D kinematics of a limited sample of 21 objects. When we 
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apply an age velocity relation to the red photometric outliers and the low gravity dwarfs we 
do find ages that are on the order of ~ 1 Gyr. We discuss the red outliers below but conclude 
that the low surface-gravity dwarfs are kinematically younger than the full or 20 pc sample. 
This result is consistent with what has already been reported through spectroscopic studies. 
There do not appear to be any low surface-gravity dwarfs flagged in the Os07 sample however 
further examination of their L and T dwarf spectra might be warranted by the discrepancy 
in ages between our samples. We suggest that kinematic studies of UCDs to date, including 
Os07, may have been plagued by small number statistics or a bias in the sample analyzed. 



6.2. Ages of the Red and Blue Outliers 

We have defined two subgroups of the ultracool dwarf population in this article that 
are both photometrically and kinematically distinct from the full or 20 pc samples. Objects 
whose J — Kg colors are sufficiently deviant are also kinematically different from the overall 
population. While we again advise caution in using the AVR, we can use it to compare the 
predicted ages of the photometric outliers to the predicted ages for the full or 20 pc samples. 
We find that the kinematics for the red outliers are consistent with a younger population 
of ultracool dwarfs whereas the kinematics for the blue outliers are consistent with an older 
population. The ~ 1 Gyr mean age for the red outliers coincides with the prediction of Os07 
for the entire L and T population. We have examined the photometry of their sample and 
concluded that the J — Kg colors of their objects are normal so the age calculation of their 
sample is not influenced by a bias from inclusion of photometric outliers. The ~ 38 Gyr 
mean age for the blue outliers is misleading. It indicates a large divergence from the full and 
20 pc samples but also indicates that the AVR must be incorrect for these objects. The more 
informative number in this case is the median T^an which, at 56 km/s, is nearly twice the 
expected value for the thin disk (see Reid & Hawley 2005). The blue photometric outliers 
most likely belong to an older population of the Galaxy such as the thick disk or the halo. 
The Wielen AVR is only valid for thin disk objects and we are unaware of an equivalent age 
relation for the halo or thick disk population. 

From our kinematic analysis we conclude that there is an age-color relation that can 
be derived for the UCD fleld population. A change in broad-band collision-induced H2 
absorption that suppresses flux at K-band is partially responsible for the near-lR color and 
consequently the age of the photometric outliers (Linsky 1969; Saumon ct al. 1994; Borysow 
et al. 1997). H2 absorption is pressure and hence gravity sensitive. Changes in H2 absorption 
effect gravity sensitive features which are used as an indicator of age. The overlap of red 
photometric outliers with low surface-gravity dwarfs and the concensus within the literature 
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that low-g dwarfs are young demonstrates the age sensitivity of H2 absorption. 

Cloud properties have also been hnked to a change in near-IR color. The analyses 
of BOS, and Gushing et al. 2008 have shown that the thickness of patchy or large-grained 
condensate clouds at the photospheres of dwarfs will lead to redder (thick clouds) or bluer 
(thin clouds) near-IR colors. The old age imphed by the kinematics of the blue outliers and 
the overlap with the UBLs suggests that there is a correlation between cloud properties and 
age or metallicity; but further investigation is warranted in this area. 

Jameson et al. (2008) have proposed a relation for inferring the ages of young L dwarfs 
using only near-infrared photometry and estimated distances. Their work supports the 
argument for an age-color relation for the ultracool dwarf population. The ages that they 
work with are no larger than ~ 0.7 Gyr. At these young ages, the surface gravities of UGDs 
change more rapidly than for ages greater than a few Gyr, so the age-color relation may be 
much stronger in the Jameson et al. (2008) sample than that seen for field dwarfs. 

7. CONCLUSIONS 

We present new proper motions for 427 late-type M, L, and T dwarfs and combine all 
previous proper motion measurements with either parallax measurements or spectrophoto- 
metric distances to compute tangential velocities for 841 M7-T9 dwarfs. We derive average 
kinematic and photometric values for individual spectral types as well as for the late-type 
M,L, and T populations as a whole. We conduct a crude U,V,W analysis and find that 
the full and 20 pc samples examined in this article have space velocities consistent with the 
Galactic thin disk population. However, there are 14 objects in the ultracool dwarf popula- 
tion that lie at the tail end of the velocity distribution and are likely to be part of an older 
Galactic population. Ages for the 20 pc sample of this kinematic study are consistent with 
the 3-6 Gyr values derived in population synthesis models; we propose that one reason for 
prior kinematic reports of ~1 Gyr mean ages for the L and T dwarf population is due to 
small number statistics or a bias in the sample analyzed. 

We find a large difference in the kinematics between red and blue photometric outliers 
and conclude that their velocity dispersions are kinematically distinct from the full or 20 pc 
samples. Analysis of the low surface-gravity and UBL subgroups also shows a distinction 
from the full and 20 pc samples. Applying an age- velocity relation we conclude that the red 
outliers and low surface-gravity subgroups are younger than the full and 20 pc samples and 
the blue outhers and UBLs are older. 
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Fig. 1. — Spectral type distribution of all late-type M, L, and T dwarfs. The overall his- 
togram is the distribution of all ultracool dwarfs in our sample. The blue shaded histogram 
shows ultracool dwarfs with proper motion measurements. The diagonally shaded histogram 
shows the distribution of ultracool dwarfs with new proper motions reported in this paper. 
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Fig. 2. — (Top): Distribution of proper motion uncertainties for the sample of 427 measure- 
ments reported in this paper. The median value is 18 mas/yr . (Bottom): Distribution of 
proper motion baselines (time between first and second epoch measurements) used in this 
survey. 
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Fig. 3. — Comparison of Right Ascension (top) and Declination (bottom) proper motion 
measured in this paper and those measured in the hterature. The straight hne represents 
a perfect agreement between measurements. The red highhghted objects are the discrepant 
proper motion measurements (see Table 4.) 
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Fig. 4. — Reduced proper motion diagram using the 2MASS J and Ks magnitudes. Late- 
type M dwarfs are marked with a black plus sign, L dwarfs are marked as a red 5 point star 
and T dwarfs are marked as blue diamonds. The line at Hk of 18 marks where M dwarfs 
are segregated from the L and T dwarfs regardless of near-IR color. This cut-off will also 
include subdwarfs and cool white dwarfs but these objects will be rare. 
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Fig. 5. — Cumulative distance distribution of all late-type M, L, and T dwarfs in our 
database. Triangles refer to the M7-M9 dwarfs, the 'X' symbols refer to all L0-L9 dwarfs 
and the plus symbols refer to all T0-T8 dwarfs. The solid line corresponds to a constant 
density distribution (N oc d^). The L and M dwarfs deviate from this distribution around 
20 pc but the T dwarfs fall off closer to 15 pc. 
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Fig. 6. — Distribution of Vtan values binned by spectral type. The top panel is the full 
astrometric sample and the bottom panel is the 20 pc sample. The asterisks refer to the 
median Vtan values and the vertical bars refer to the standard deviation or dispersion of 
velocities. 
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Fig. 7. — The overall histogram is the tangential velocity distribution for the entire sample 
and the diagonally shaded histogram is the 20 pc sample. Both Vtan distributions peak in 
the 10-30 km s~^ bins. 
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Fig. 8. — Histogram of U,V,W velocities. Plotted for each velocity is 1) each object in the 
astrometric sample (large histogram) 2) a 30 degree restriction on objects and 3) a 60 degree 
restriction (smallest histogram). The 30 and 60 degree restrictions are placed on the X,Y or 
Z axis and correspond to removing the U,V, or W velocity respectively for objects in cones 
of noted radius around the respective axis. 
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Fig. 9. — J — Kg colors of late-type dwarfs. We compute the average values for each spectral 
type (binned by 1 subtype) from the 2MASS photometry of a select sample of dwarfs and 
then flag objects as photometric outliers when they are either twice the standard deviation 
oi J — Kg or 0.4 magnitude redder or bluer than the average value. Red symbols above the 
plotted range of J — Kg colors are the red outliers and blue symbols below the plotted range 
oi J — Kg colors are the blue outliers. 
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Fig. 10. — The spread of tangential velocities for objects marked as red outliers (top panel) 
and blue outliers (bottom panel). The red population has a fairly tight dispersion and the 
blue population has a fairly wide dispersion compared to the full sample suggesting a link 
between near-IR color and age. The dashed line in each plot represents the median Vtan 
value for the outlier group and the solid black hues represent the dispersion. 
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Fig. 11. — A scatter plot showing Vtan as a function of the deviation in J — Kg color from the 
average at a given spectral type. The blue outhers appear to move faster on average than 
the red outliers. To demonstrate this we have over-plotted the average Vtan with dispersion 
for the blue and red photometric outliers as well as for the full astrometric sample (dashed 
lines). 
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Fig. 12. — Top: A plot of median Vtot and atot values calculated from the U,V,W velocities 
for the 20 pc sample of F through T objects. Bottom: A plot of median Vtan and atan 
values calculated from the proper motions and distances for the 20 pc sample of F through 
T objects. 
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Table 1 . Properties of Instruments Used for Astrometric Measurements 



Telescope 


Instrument 


Band 


FOV 
(arcminutes) 


Plate Scale 
(arcsecond / pixel) 


Dates 


Seeing 
(arcseconds) 


Sources Observed 


CTIO 0.9m 


CFIM 


I 


4.5 


0.40 


2007 Sep 23 - 26 


0.8 - 2.0 


42 


MDM 1.3m 


TIFKAM 


J 


5.6 


0.55 


2007 Nov 20 - 24 


1.0 - 2.5 


66 












2008 Apr 22 - 28 


0.8 - 2.5 


80 


CTIO 1.5m 


CPAPIR 


J 


35.0 


1.02 


2005 Oct 19 


1.0 - 1.5 


4 












2006 Aug 21 


1.0 - 2.0 


7 












2007 Mar 04 


0.9 - 1.8 


28 












2007 Mar 23 


0.9 - 1.3 


39 












2007 Dec 03 - 06 


1.0 - 2.0 


35 












2008 Jan 15 - 23 


0.7 - 2.5 


248 
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Table 4. Discrepant Proper Motion Values 



Name 


fJ,aCOS(0) 




^ciCos(d) 


^dec 


Reference 




C'/yr) 


C'/yr) 


C'/yr) 


C'/yr) 






This Paper 


This Paper 


Literature 


Literature 




SIPS J0050-1538 


-0.229 ± 0.018 


-0.494 ± 0.019 


-0.495 ± 0.039 


-0.457 ± 0.038 


16 


2MASSJ02271036-1624479 


0.426 ± 0.016 


-0.297 ± 0.017 


0.509 ± 0.016 


-0.303 ± 0.010 


16 


2MASSJ09393548-2448279 


0.592 ± 0.019 


-1.064 ± 0.021 


0.486 ± 0.031 


-1.042 ± 0.055 


41 


zivi^oo vv u J. xuuoyo o i ^ i oo 


n n^n -}- n ni 9 




U.-LXO Zn \J .KJKJO 


-H j<6i -1- n n^Q 

u.ouj- zn u.uof 


16 








0.013 ± 0.015 


-0.778 ± 0.013 


9 








0.06 ± 0.04 


-0.82 ± 0.07 


36 


2MASSJ13411160-3052505 


0.030 ± 0.013 


-0.134 ± 0.015 


0.109 ± 0.014 


-0.163 ± 0.022 


17 


2MASSJ13475911- 7610054 


0.203 ± 0.005 


0.038 ± 0.020 


0.257 ± 0.063 


0.287 ± 0.063 


22 








0.193 ± 0.011 


0.049 ± 0.019 


35 


2MASSWJ1448256+103159 


0.262 ± 0.022 


-0.120 ± 0.022 


0.70 ± 0.15 


-0.10 ± 0.16 


36 








0.249 ± 0.015 


-0.099 ± 0.016 


10 


2MASSWJ1507476-162738 


-0.128 ± 0.014 


-0.906 ± 0.015 


-0.043 ± 0.011 


-1.037 ± 0.255 


16 








-0.1615 ± 0.0016 


-0.8885 ± 0.0006 


15 








-0.147 ± 0.003 


-0.890 ± 0.002 


12 








-0.09 ± 0.11 


-0.88 ± 0.06 


36 


2MASSJ15485834-1636018 


-0.210 ± 0.016 


-0.107 ± 0.017 


-0.189 ± 0.016 


-0.176 ± 0.015 


17 








-0.098 ± 0.043 


-0.161 ± 0.042 


22 


2MASSWJ1555157-095605 


0.950 ± 0.015 


-0.767 ± 0.015 


0.929 ± 0.014 


-2.376 ± 0.017 


10 








0.961 ± 0.017 


-0.835 ± 0.014 


16 








-0.400 ± 1.200 


-1.900 ± 1.100 


9 


2MASSJ19360187-5502322 


0.169 ± 0.009 


-0.298 ± 0.016 


0.603 ± 0.037 


-0.579 ± 0.035 


16 








0.22 ± 0.29 


-0.19 ± 0.28 


36 


2MASSJ22551861-5713056 


-0.216 ± 0.011 


-0.260 ± 0.020 


0.394 ± 0.321 


-1.525 ± 0.319 


22 








-0.16 ± 0.11 


-0.32 ± 0.13 


36 


2MASSJ23302258-0347189 


0.223 ± 0.022 


0.014 ± 0.022 


0.349 ± 0.051 


-0.107 ± 0.016 


16 








0.232 ± 0.017 


0.032 ± 0.013 


10 





Note. — Details on the discrepant proper motion objects. We note only objects whose proper motion values were 
discrepant by more than 2cr. Proper motion references are listed in Table 3. 



-77- 



Table 5. Median Photometric and Kinematic Properties of Ultracool Dwarfs 



SpT 








Median Dist 










NRed 








{"M 


C'/yr) 


(pc) 


(pc) 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


M7 


88 


0.261 


0.553 


25 


10 


160 


1.08 


0.19 





1 


M8 


114 


0.210 


0.403 


23 


8 


147 


1.14 


0.18 


1 


1 


M9 


71 


0.204 


0.357 


22 


10 


107 


1.20 


0.22 


1 





LO 


93 


0.111 


0.211 


32 


19 


92 


1.31 


0.37 


4 


1 


LI 


83 


0.208 


0.301 


31 


21 


82 


1.39 


0.37 


4 


1 


L2 


58 


0.185 


0.209 


32 


17 


63 


1.52 


0.40 


5 


1 


L3 


64 


0.189 


0.398 


33 


17 


67 


1.65 


0.39 


1 


1 


L4 


50 


0.183 


0.284 


27 


12 


44 


1.73 


0.40 


2 


2 


L5 


43 


0.323 


0.281 


24 


12 


43 


1.74 


0.40 





1 


L6 


36 


0.215 


0.339 


26 


12 


31 


1.75 


0.40 


4 


2 


L7 


21 


0.247 


0.186 


23 


9 


15 


1.81 


0.40 





2 


L8 


16 


0.280 


0.368 


19 


8 


16 


1.77 


0.33 


2 





L9 


3 


0.424 


0.200 


20 


6 


7 


1.69 


0.19 








TO 


9 


0.333 


0.165 


18 


4 


8 


1.63 


0.40 








Tl 


11 


0.289 


1.336 


23 


9 


10 


1.31 


0.40 


1 


1 


T2 


13 


0.350 


0.285 


15 


7 


15 


1.02 


0.40 


1 





T3 


7 


0.183 


0.135 


26 


6 


5 


0.63 


0.40 


1 





T4 


13 


0.323 


0.219 


23 


9 


6 


0.26 


0.40 








T5 


20 


0.340 


0.351 


15 


3 


12 


0.07 


0.39 








T6 


15 


0.594 


1.217 


11 


18 


5 


-0.30 


0.40 


2 


1 


T7-T8 


13 


1.218 


0.764 


9 


3 


6 


-0.08 


0.40 





1 


M7-M9 


273 


0.222 


0.445 


23 


9 


414 


1.12 


0.22 


2 


2 


L0-L9 


467 


0.189 


0.292 


29 


17 


460 


1.53 


0.40 


22 


11 


T0-T9 


101 


0.373 


0.801 


15 


10 


67 


0.74 


0.40 


5 


3 





Note. — To calculate the {J — Ks)avg for each spectral type, we chose only objects that were not identified as binaries,young 
cluster members, subdwarfs and/or had aj and < 0.20. 
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Table 6. 20 pc Sample 



SpT 


N 


N 

High Vtan 


median Vtan 


median Vtan 
with high Vtan 


Ctan 


0"tan 

with high Vtan 


Age 


Age 
with high Vtan 








(km s-i) 


(km s-l) 


(km s-1) 


(km s-l) 


(Gyr) 


(Gyr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


M7 


29 





25 


25 


20 


20 






M8 


37 


1 


33 


33 


20 


25 






iviy 


z / 


1 
1 


ZD 




oo 
ZZ 


ZO 






LO 


9 





19 


19 


21 


21 






LI 


19 





30 


30 


29 


29 






L2 


10 





27 


27 


16 


16 






L3 


12 


3 


32 


38 


20 


46 






L4 


15 


1 


27 


27 


20 


28 






L5 


16 





27 


27 


21 


21 






L6 


10 





28 


28 


24 


24 






L7 


9 





30 


30 


9 


9 






L8 


12 





25 


25 


20 


20 






L9 


2 





41 


41 












TO 


6 





32 


32 


15 


15 






Tl 


3 





66 


66 


28 


28 






T2 


8 





26 


26 


5 


5 






T3 


1 





39 


39 












T4 


5 





21 


21 


16 


16 






T5 


20 





21 


21 


23 


23 






T6 


14 





44 


44 


22 


22 






T7 


10 


1 


45 


54 


15 


34 






T8 


3 





57 


57 


8 


8 






M7-M9 


93 


2 


29 


29 


21 


24 




n+i.7 
J-U_i 4 


L0-L9 


114 


5 


27 


27 


21 


26 


o 2+1-1 




T0-T9 


70 


1 


30 


31 


20 


24 


^•°-0.8 


4 6+1-6 

^•"^-1.3 



Note. — The age range is calculated from the Wielen (1977) age-velocity relation for the disk which uses a value of a of 
(1/3). 
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Table 7. Full Astrometric Sample 



SpT 


N 


N 

High Vtan 


median Vtan 


median Vtan 
with high Vtan 


Ctan 


0"tan 

with high Vtan 


Age 


Age 
with high Vtan 








(km s-i) 


(km s-l) 


(km s-1) 


(km s-l) 


(Gyr) 


(Gyr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


M7 


88 





27 


27 


19 


19 






M8 


114 


1 


27 


27 


21 


23 






iviy 




1 
1 


ZO 


Zo 


1 n 

iy 


ZL 






LO 


93 


1 


19 


19 


16 


21 






LI 


83 


2 


32 


33 


23 


27 






L2 


58 





26 


26 


18 


18 






L3 


64 


3 


30 


32 


18 


27 






L4 


50 


1 


25 


27 


20 


23 






L5 


43 





25 


25 


20 


20 






L6 


36 


1 


26 


27 


18 


24 






L7 


21 


1 


28 


28 


13 


22 






L8 


16 





25 


25 


19 


19 






L9 


3 





38 


38 


17 


17 






TO 


9 





26 


26 


13 


13 






Tl 


11 





31 


31 


25 


25 






T2 


13 





26 


26 


11 


11 






T3 


7 





25 


25 


10 


10 






T4 


13 





32 


32 


22 


22 






T5 


20 





21 


21 


23 


23 






T6 


15 





36 


36 


23 


23 






T7 


10 


1 


45 


54 


15 


34 






T8 


3 





57 


57 


8 


8 






M7-M9 


273 


3 


26 


26 


19 


21 


9 c+lt.H 


■''•^-0.9 


L0-L9 


467 


10 


26 


26 


19 


23 


9 t:+0.9 




T0-T9 


101 


1 


29 


29 


20 


23 


9 7+1.0 
^•'-0.8 


4.ot;:t 



Note. — The age range is calculated from the Wielen (1977) age-velocity relation for the disk which uses a value of a of 
(1/3). 
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Table 8. Average Kinematics and Ages for the Subgroups 



SpT N Median Vt an (^tan -Agc Range 

km km Gyr 

(1) (2) (3) (4) (5) 



M7-T9/BLUE 


16 


53 


47 


Q7 q+12-6 


M7-T9/RED 


29 


26 


16 


1 n+0.5 

-■-■^-0.4 


M7-L9/BLUE 


13 


56 


50 


46.0^;^-^ 


M7-L9/RED 


24 


26 


15 


1 n+°* 


UBLs 


10 


99 


47 


•^'•^-10. 3 


Low Gravity 


26 


18 


15 


1 o+'^-'' 





Note. — Tlie age range is calculated from the Wielen (1977) age- 
velocity relation for the disk which uses a value of a of (1/3). 



Table 9. Details on Red Photometric Outliers 



Source Name 2MASS J 2MASS Ks HqCos((5) jjj /a Ref. SpT SpT Vtan Note' 

(mag) (mag) ("/yr) ("/yr) (opt) (IR) (kms"!) 

(1) (2) (3) (4) (6) (6) (7) (8) (9) (11) 



2MASS ,700374306-5846229 


15. 


37 


± 





.05 


13, 


,59 


± 


0, 


,05 


0.049 ± 0.010 


-0.051 ± 0.020 


19 


LO 






18 


± 5 


LG 


SDSSp J010752. 33-1-004156.1 


15. 


,82 


± 





.06 


13, 


,71 


± 


0, 


.04 


0.628 ± 0.007 


0.091 ± 0.004 


44 


L8 




L5.5 


46. 


,9 ± 3.3 




2MASS J01244B99-S746379 


16. 


,31 


± 


0, 


.11 


14, 


.32 


± 


0, 


.09 


-0.003 ± 0.010 


0.018 ± 0.019 


19 


LO 






7 ± 7 


LG 


2MASS J01415823-4633574 


14. 


,83 


± 





.04 


13, 


,10 


± 


0, 


.03 


0.104 ± 0.017 


-0.026 ± 0.024 


19 


LO 




LO 


21 


± 4 


LG 


2MASS ,701490895+2956131 


13, 


,45 


± 





.02 


11 


,98 


± 


0, 


,02 


0.1757 ± 0.0008 


-0.4021 ± 0.0007 


15 


M9 


.5 




46. 


8 ± 0.7 




2MASSI ,70243137-245329 


15. 


,42 


± 





.06'= 


15 


,22 


± 





,06" 


-0.288 ± 0.004 


-0.208 ± 0.003 


44 






T6 


IS. 


± 0.7 




2MASS J03231002-4631237 


15. 


39 


± 





.07 


13, 


,70 


± 


0, 


.05 


0.060 ± 0.013 


-0.010 ± 0.019 


19 


LO 






16 


± 4 


LG 


2MASS J03264225-2102()57 


16. 


13 


± 





.09 


13 


,92 


± 


0, 


.07 


0.108 ± 0.014 


-0.146 ± 0.015 


19 


L4 






36 


± 6 




2MASS J03421621-6817321 


16. 


,85 


± 





.14 


14, 


,54 


± 


0, 


.09 


0.064 ± 0.007 


0.021 ± 0.018 


19 


L2 






26 


± 6 




2MASS J03552337-I-1133437 


14. 


,05 


± 


0, 


.02 


11, 


.53 


± 


0, 


.02 


0.192 ± 0.017 


-0.613 ± 0.017 


19 


L6 






26 


± 6 


LG 


2IVIASS J0436 1456- 1414468 


11. 


,88 


± 





.03 


9.95 ± 0.02 


0.009 ± 0.014 


0.016 ± 0.014 


19 


M8 






1 ± 1 


LG 


2MASS J05012406-0010452 


14, 


,98 


± 





.04 


12 


96 


± 





,04 


0.158 ± 0.014 


-0.139 ± 0.014 


19 


L4 






24 


± 5 


LG 


2MASSI ,70512063-294954 


15. 


,46 


± 





.06 


13 


29 


± 





,04 


-0.028 ± 0.016 


0.099 ± 0.018 


19 


L4. 


6 




13 


± 3 




2MASS ,705361998-1920396 


15, 


77 


± 





.08 


13 


,85 


± 





,06 


0.017 ± 0.017 


-0.024 ± 0.018 


19 


LI 






8 ± 5 




AB Pic b 


16, 


18 


± 





.10 


14, 


,14 


± 





,08 


0.0141 ± 0.0008 


0.0452 ± 0.0010 


34 






LI 


10. 


2 ± 0.4 


VLMC 


SDSS J080959. 01-1-443422. 2 


16. 


,51 


± 





.06" 


14, 


,34 


± 


0, 


.06" 


-0.198 ± 0.014 


-0.214 ± 0.019 


19 






L6 


35 


± 7 




SDSS J085834.42-H32B627.7 


16. 


,62 


± 


0, 


.06* 


14, 


.69 


± 


0, 


.06" 


-0.760 ± 0.023 


0.075 ± 0.023 


19 






Tl 


66 


± 3 




G 196-3B 


14. 


,83 


± 





.05 


12, 


,78 


± 


0, 


.03 


-0.133 ± 0.040 


-0.185 ± 0.015 


10 


L2 






36 


± 6 


VLMC 


2MASS J12123389+0206280 


16, 


13 


± 





.13 


14, 


,19 


± 


0, 


,09 


0.065 ± 0.021 


-0.141 ± 0.021 


19 






LI 


49 


± 9 




2MASS ,713243559+6358284 


15, 


,60 


± 





.07 


14, 


,06 


± 





,06 


-0.343 ± 0.064 


-0.260 ± 0.048 


32 






T2 


26 


± 6 




SDSSp .7132629.82-003831.5 


16, 


,37 


± 





.06" 


14, 


,17 


± 





,06" 


-0.226 ± 0.008 


-0.107 ± 0.006 


44 


L8 




L5.5 


23. 


8 ± 3.2 




SDSS J141530. 05 + 572428. 7 


16, 


72 


± 





.06" 


15 


,49 


± 





,06" 


0.043 ± 0.013 


-0.345 ± 0.025 


19 






T3 


36 


± 12 




2MASS J15311344+1641282 


15. 


,58 


± 





.06 


13, 


,80 


± 


0, 


.05 


-0.076 ± 0.025 


0.040 ± 0.026 


19 






LI 


21 


± 8 




2IV1ASSI J1726000+153819 


16. 


,67 


± 


0, 


.07 


13, 


.66 


± 


0, 


.06 


-0.031 ± 0.013 


-0.048 ± 0.014 


10 


L2 






13 


± 3 




SDSS J175805. 46+463311. 9 


16. 


17 


± 


0, 


.06= 


16, 


.99 


± 


0, 


.06= 


0.026 ± 0.016 


0.594 ± 0.016 


10 






T6.6 


34 


± 4 




2MASS J21481633+4003594 


14, 


15 


± 





.03 


11 


,77 


± 


0, 


.02 


0.770 ± 0.018 


0.456 ± 0.024 


19 


L6. 


6 




30 


± 6 




2MASS J21512543-2441000 


15, 


,75 


± 





.08 


13, 


,65 


± 


0, 


.05 


0.278 ± 0.014 


-0.021 ± 0.015 


19 


L3 






56 


± 6 




2MASSW J2206450-421721 


15. 


,56 


± 





.07 


13, 


,61 


± 


0, 


.06 


0.111 ± 0.013 


-0.182 ± 0.018 


19 


L2 






46 


± 6 




2MASSW J2244316+204343 


16. 


,47 


± 


0, 


.06= 


13, 


.93 


± 


0, 


.06" 


0.252 ± 0.014 


-0.214 ± 0.011 


10 


L6. 


6 


L7.B 


30 


± 3 





Note. — See Table 3 for references and notes referred to in this table. 
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Table 10. Details on Blue Photometric Outliers 



Source Name 


2MASS J 


2MASS K, 


/ic(Cos(<5) 




^ Ref. 


SpT 


SpT 


Vton 


Note' 




(mag) 


(mag 


) 


C'/yr) 


CM) 




(opt) 


(IR) 


(km s-1) 




(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(11) 


HD 3661B 


16.16 ± 0.03 


16.87 ± 0.06 


-0.4611 ± 0.0007 


-0.3709 ± 0.0007 


34 




T7.6 


31.2 ± 0.3 


VLMC 


SSSPM J0134-6316 


14.61 ± 0.04 


13.70 ± 0.04 


0.077 ± 0.008 


-0.081 ± 0.009 


30 




LO 


19 ± 2 




2MASS J02530084+1652532 


8.39 ± 0.03 


7.59 ± 0.05 


3.404 ± 0.005 


-3.807 ± 0.005 


23 


M7 




92.9 ± 1.0 




SDSS ,7090900.73+652527.2 


16.00 ± 0.06" 


15. IB ± 


O.OB" 


-0.217 ± 0.003 


-0.138 ± 0.008 


19 




Tl 


28 ± 1 




2MASS J09211410-2104446 


12.78 ± 0.02 


ll.fiH ± 


0.02 


0.244 ± O.OIB 


-0.908 ± 0.017 


19 


L2 




56 ± 4 


UBL 


SDSS J093109. 56+032732. 5 


16.75 ± O.IO'' 


15.fl5 ± 


o.io'- 


-0.612 ± 0.018 


-0.131 ± 0.018 


19 




L7.5 


108 ± 23 


UBL 


2MASSI J0937347+293142 


14.58 ± 0.06^= 


15.51 ± 


0.12'^ 


0.973 ± 0.005 


-1.298 ± 0.006 


44 


d/sdT6 


T6 


47.2 ± 1.1 




SDSS J103321.92+400549.6 


16.88 ± 0.06'* 


16.63 ± 


0.10" 


0.154 ± 0.013 


-0.188 ± 0.018 


19 




L6 


53 ± 10 


UBL 


SDSS J112118. 67+433246. B 


17.19 ± 0.10* 


16.16 ± 


0.08" 


-0.067 ± 0.024 


0.026 ± 0.033 


19 




L7.5 


14 ± 6 


UBL 


2MASS J11263991-5003550 


14.00 ± 0.03 


12.83 ± 


0.03 


-1.570 ± 0.004 


0.438 ± 0.011 


20 


L4.5 


L9 


lOB ± 11 


UBL 


SDSS J114805. 02+020350. 9 


15.52 ± 0.07 


14.51 ± 


0.12 


0.237 ± 0.026 


-0.322 ± 0.013 


10 


LI 




96 ± 8 




2MASS J12162161+4456340 


16.35 ± 0.10 


15.02 ± 


0.12 


-0.035 ± 0.014 


-0.004 ± 0.019 


19 


L5 




7 ± 3 




SDSS J142227.25+221557.1 


17.01 ± 0.06*' 


16.67 ± 


0.06" 


0.047 ± 0.019 


-0.054 ± 0.020 


19 




L6.5 


14 ± 6 


UBL 


DENIS-P J170548. 38-051646. 7 


13.31 ± 0.03 


12.03 ± 


0.02 


0.129 ± 0.014 


-0.103 ± 0.016 


10 




L4 


9 ± 1 




2MASSI J1721039+334416 


13.63 ± 0.02 


12.49 ± 


0.02 


-1.864 ± 0.017 


0.602 ± 0.017 


10 


L3 




144 ± 13 


UBL 


2MASS J18261131+3014201 


11.66 ± 0.02 


10.81 ± 


0.02 


-2.280 ± 0.010 


-0.684 ± 0.010 


28 


M8.6 




132 ± 9 





Note. — 



See Table 3 for references and notes referred to in this table. 
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